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Exact analytic solutions are found for the Dirac equation in 2+1 dimensions for a spin-one- 
half particle in a combination of the Lorentz 3-vector and scalar Coulomb as well as Aharonov- 
Bohm potentials. We employ the two-component Dirac equation which contains a new parameter 
introduced by Hagen to describe the spin of the spin-1/2 particle. We derive a transcendental 
equations that implicitly determine the energy spectrum of an electron near the negative-energy 
continuum boundary and the critical charges for some electron states. Fermion pair production 
from a vacuum by a strong Coulomb field in the presence of the magnetic flux tube of zero radius is 
considered. It is shown that the presence of the Ahanorov-Bohm flux tends to stabilize the system. 
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I. INTRODUCTION 

The Aharonov-Bohm (AB) effect is one of the most intriguing effects of a truly quantal nature Ever since its 
discovery, the AB effect has been analyzed in various physical situations in numerous works [|J . In recent years there 
has been considerable interest in the problem of scattering of spin-1/2 particle off an AB potential in 2+1 dimensions. 
The results of [l[ for nonrelativistic case modified by using the Dirac equation in 2+1 dimensions were applied to 
many problems. For instance, solutions to the two-component Dirac equation in the AB potential were first discussed 
by Afford and Wilczek in Q in a study of the interaction of cosmic strings with matter. In particular^ a mechanism 
of particle production due to the nonstatic AB potential of the moving cosmic string was discussed in [3| . Relativistic 
quantum AB effect was studied in Ref. [J] for the free and bound fermion states by means of exact analytic solutions 
of the Dirac equation in 2+1 dimensions for a combination of AB potential, Lorentz three- vector and scalar Coulomb 
potentials. 

In [5[ the scattering of spin-polarized fermions in an AB potential was considered in 2+1 dimensions. There the 
particle spin was introduced into the two-component Dirac equation as a new parameter. The term including this 
new parameter appears in the form of an additional delta-function interaction of spin with magnetic field in the Dirac 
equation. Solutions of the Dirac equation were then interpreted for the case of 3+1 dimensions. Similar problems were 
also discussed in [6j by taking the AB flux tube with a small but finite radius in (3+l)-dimensional Pauli equation, 
and in (2+l)-dimensional Dirac equation. 

In this paper we would like to study how various physical fields affect the properties of a bound Dirac fermion in 
2+1 dimensions. Specifically, we study the energy spectrum of the fermion in a combination of the Lorentz 3-vector 
and scalar Coulomb as well as AB potentials. We also consider the influence of the magnetic flux tube of zero radius 
on the so-called critical charge which determines the onset of instability of the system, and on the probability of 
production of an electron-positron pair from a vacuum by a strong vector Coulomb field. We note here that in 3+1 
dimensions, analytic solutions of such problems, even for the Schrodinger equation in the Coulomb and AB potentials, 
have not yet been found. 

This paper is organized as follows. In section II we find the exact bound states solutions of the Dirac equation in 
2+1 dimensions for a combination of the Lorentz 3-vector, scalar Coulomb and AB potentials for the spin-1/2 particle. 
The formulas for eigen-energies of the relativistic fermion are obtained and discussed. In section III the critical charge 
and the probability of production of an electron-positron pair from a vacuum by a strong vector Coulomb field are 
calculated. Using a simplified model with a truncated Coulomb potential we show that the critical charge and the 
production probability of pair are influenced by the magnetic flux and spin particle. 

II. ENERGY SPECTRUM OF A DIRAC FERMION IN PLANAR COULOMB AND 

AHARONOV-BOHM POTENTIALS 

The (2+l)-dimensional Dirac equation of a fermion of mass m and charge e = — eo < in a vector potential 
and a Lorentz scalar potential U is (c = K = 1) 



(7^-m-C/)* = 0, 



(1) 
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where P^ = —id^ — eA^ is the generalized fermion momentum operator. The Dirac 7 M matrices are conveniently 
defined in terms of the Pauli spin matrices as 

7° = 0-3, 7 1 = isai, 7 2 = ia 2 , (2) 

Following [f|, here s is a new parameter characterizing twice the spin value s = ±1 for spin "up" and "down", 
respectively. 

We are interested in finding exact analytic solutions of the Dirac equation for both signs of s in an AB potential, 
which is specified in the Cartesian or cylindrical coordinates as 

A = 0, A x — j", Ay — — J - ; A" = 0, A r = 0, A^ = — , 



r = \/ x 2 + y 2 , ip = tan 1 (y/x) (3) 
and a Lorentz 3-vector potential (AZ'(r)) and a scalar potential(J7(r)) potential defined by 



A°(r) = -^-, A r = 0, A tp = 0; U(r) = -- ( O> 6>0) (4) 
\e\r ' r 



In [1] only the case for s — 1 was considered. Note that if e -B = N, where N is an integer, then the magnetic field 
flux is quantized as $ = where $0 = 27r/eo is the elementary magnetic flux called the "fluxon". 

The Dirac Hamiltonian for this system is 

H D = aiP 2 - sa 2 P 1 +o 3 (m + U{r)) - e A a (r). (5) 

The vector potential A^ in the generalized momentum P^ is the sum of the AB potential ([3]) and the Lorentz 3-vector 
potential ([3]). The total angular momentum J z = L z + sa^,/2, where L z = —id/ dtp, is a conserved quantity. 
We seek positive energy solutions of Eq. fl} in the form (see, also, Ref. @, HQ) 

$(t, x) = exp(-iEt + il(p)ip(r, <p) , (6) 
V 2tt 



where E > is the fermion energy, I is an integer, and if)(r, ip) is a two-component function (i.e. & 2-spinor) 

mw) = ( JJzL)- ( 7 ) 



?(r)e 1 ' 

The wave function ^ is an eigenfunction of the conserved total angular momentum J z with eigenvalue j — I + s/2. 
From the Dirac equation one finds that f(r) and g(r) satisfy the following equations : 

s- / + IE + m -\ g = 0, 

ar r \ r ) 

dg l + s + e B ( a + b\ 
s— H g-[E-m-\ 7=0. (8) 



dr 



r 



Compared with the situation studied in H in which only the vector Coulomb potential (a) is present, one notes that 
the effect of the AB potential (B) appears only in modifying the angular momentum I. As in [8(, we now assume /(r) 
and g(r) to have the form (see, for example, (IP ]) 

f{r) = V^+Ee- x / 2 x^(Qt + Q s 2 ) , 

g{r) = ^~Ee- x ' 2 x^ (Q\ - Q|) , (9) 

where 



x = 2Ar, A = y/m 2 - E 2 , (10) 

and 



l -±\J(l + e B+ S -) 2 -a 2 + b 2 (11) 
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determines the asymptotic behavior of the wave function for small r. 

From Eq. (fTl"]) one sees that when a 2 < (I + e B + s/2) 2 + b 2 the quantity 7 S is real, and must be chosen positive 
to ensure normalizability of the wave function. If a 2 > (I + cqB + s/2) 2 + b 2 then the two roots of 7s are imagi nary 
and the corresponding wave functions oscillate as r — > 0, which indicates the occurrence of Klein's paradox [10L 
We shall consider this situation in the next section. In what follows we shall take 



7a 



(l + e B 



b 2 . 



(12) 



For wave functions that are finite at x = 0, the functions Q\ 2 for s = ±1 are given by the confluent hypergeometric 
function F(a±, C\;x): 



Ql = AF I 7s 
Q s 2 = CF ( 7s 
The constants A and C are related by 



s aE + mb 



C 



(s7 s + s/2) - (Ea + mb) /A 
I + e B + s/2 + (ma + bE)/X 



- 2;x 
1-2; x 

A. 



(13) 



(14) 



The wave function is normalizable if both the hypergeometric functions Q\ and Q| are reduced to polynomials. 
For this the parameter a of F(a, c; x) must be a negative integer or zero. Denoting 



Is 



s 

l ~2 



Ea + mb 



(15) 



one can see that if n r — 1,2,3,.. ., then Q\ and Q\ are reduced to polynomials for s = 1. If n r = 0, then only Q\ is 
reduced to a polynomial. But the relation n r = implies that 



1 Ea + mb 



7i+2 



A 



> 0. 



(16) 



Then it follows from (Tl6l) and the relation 



7s + 2 ' 



A 2 



i 

-E 2 = (z + eo£+! 



that 



Eb + ma 
A 



l + e B- 



(17) 



(18) 



If / + eoB + 1/2 > 0, then C = 0, hence Q\ = and the required condition is not violated (Q\ is a polynomial). If 
I + eoB + 1/2 < 0, then C = A, and Q\ is still a divergent function. The following values of n r are hence admissible: 
0,1,2,... for l + e B + l/2 > 0, and 1,2,... for l + e B + l/2 < 0. 

For s = — 1, Q^ 1 and Q^ 1 are reduced to polynomials if n r — 0, 1, 2, 3, . . .. But in this case n r also can be equal to 
— 1. If n r = —1, then only Q^ 1 is reduced to a polynomial. But for n r = —1 from (|15p one obtains 



and Eq. (jTSJ) becomes 



7-1 + 2 



Eb + ma 



A 



Ea + mb 
A 



I + e B 



> 



It is convenient to rewrite the relation (fT4|) for s = — 1 as 

A _ I + e B- 1/2 + (ma + bE)/X 
(7_i + 1/2) + (Ea + mb)/X 



C. 



(19) 



(20) 



(21) 
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So, if I + cqB — 1/2 < 0, then A — 0, hence Q 1 1 = and the required condition is not violated (Q 2 1 is a polynomial). 
If I + e B — 1/2 > 0, then A = —C, and Q^ 1 is divergent. Hence the following values of n r are admissible for s = — 1: 
-1,0,1,2,... for Z + e a B- 1/2 < 0, and 0,1,2,.. . for I + e Q B - 1/2 > 0. 
We can rewrite the equation (|15[) for the energy spectrum as 



= n r + 7s 



s 

l -2 



u. 



from which we obtain finally the discrete fermion energy levels in the form 



ab 




(22) 



(23) 



One sees that the energy spectrum is influenced by the magnetic flux through j s given by Eq. (|12p. On the other 
hand, for the flux that is integer in the unit $o the energy spectrum is the same as in the absence of magnetic flux. 
The spectrum of the system changes only for flux that is not integer in the unit $o. For such magnetic fluxes the 
energy spectrum is likely to be observed just as the AB effect. For flux that is not integer or half integer all the energy 



levels are doubly degenerate; the levels with Z, 



1. 



T and I + 1, n r , s 



-1 coincide. This reflects the fact 



that the fermion energy does not depend upon spin in the field configuration considered. 
If the scalar Coulomb potential is absent, then the energy spectrum is given by 



E n i = m 



(TV + 1/2 - s/2 + y/(l + e B + s/2) 2 - a 2 ) 2 



-1/2 



(24) 



This expression makes sense only when \l + e$B + s/2\ > a, a condition that forbids the existence of the energy levels 
with I + e B + s/2 = 0. 

In the nonrelativistic Schrodinger limit, the expression for the energy spectrum becomes 



E 



2{n r + 1/2 - s/2 +\l + e B + s/2\) 2 



(25) 



For flux eoB that is not integer or half integer it has a full analogy with the Rydberg correction. Note that for flux 
that is half integer in the unit $o the energy spectrum (1251) depends only on integer number. It can be observed using 
spectroscopy. It is interesting that for such magnetic fluxes the cross section in the AB scattering is maximal. It is 
seen that the eigen-energies of a fermion in these electromagnetic field combination are periodic function of the flux 
like the case of the motion of fermion in a closed ring in the absence of the two-dimensional Coulomb potential. The 
energy spectrum (f2"5j) repeats itself every time when the change in the flux eoB is integral. 



III. CRITICAL CHARGE AND PAIR PRODUCTION IN VECTOR COULOMB AND 

AHARONOV-BOHM POTENTIAL 

We now consider the problem of stability of the system considered in the previous section when a becomes large. 
For simplicity we ignore the scalar Coulomb potential in this section, i.e. we set 6 = 0. Such consideration is relevant 
to the stability of the vacuum of quantum electrodynamics in a strong Coulomb field. In the absence also of the 
AB potential this problem in 3 + 1 dimensions had been extensively studied in [l2|, [l3|, [TJ, 0j| [lj| 03, [3- The 
corresponding system in 2 + 1 dimensions was considered in @, Q . Now we would like to see how the presence of the 
AB potential may affect the stability of the system. 

From Eq. (|2"3"]) we see that the lowest electron states in the vector Coulomb and AB potentials are those with 
n r = I = 0, s = 1 and n r = — 1, / = 0, s = — 1. The electron energy in these states can be written as 



{e,B + s/2f-a? 

Es ~ m i (e B + s/2r ■ (26) 

For definiteness, in this section we shall consider positive flux (e^B > 0). The case with negative flux can be discussed 
similarly with the signs of I and s flipped: it is just the mirror image of the case with positive flux with respective to the 
xy-plane. For magnetic flux which is positive and half integer in $0; the energy of the state with n r = — 1, 1 = Q, s = — 1 
is divergent and imaginary for any value of a. So it is reliable to assume that the electron ground state in this case is 
the state with n r — I — and s = 1 . 



5 



The energy Eo(a) of this state as a function of a becomes zero at a = a^{B) = eoB + 1/2 and purely imaginary 
when a > a cr (B). This implies that E Q (a) becomes meaningless at a > ao(B), and the wave function oscillates with 
infinite frequency at the origin as mentioned in Sect. 2. It may appear that one cannot determine the spectrum 
beyond ao(B). However, one notes that in reality the source of the Coulomb field has finite extension so that the 
potential remains finite at the origin. In this case the wave function is regular at the origin and the energy levels can 
be traced continuously beyond ao(B). 

The problem of a (2 + l)-dimensional Dirac particle in a strong field of a truncated (at small distances) Coulomb 
potential in the absence of AB potential was considered in Q, and the expression for the electron energy spectrum 
was obtained in [9(. There it was shown that as the strength of the Coulomb source a increases, the lowest energy 
level (c.f. Eq. @) E S=1 {B = 0) = m\/l - 4a 2 was pulled towards the negative continuum. This energy becomes 
negative for a 2 > 1/4 (in the truncated potential) and may reach the negative-energy continuum boundary —m. 
When E s —x(B = 0) dives into the negative continuum the vacuum of quantum electrodynamics becomes unstable and 
particle-antiparticle pair is created spontaneously. The value a — a cr for which the lowest energy level coincides with 
— m is called the critical charge for the ground state. 

Let us consider the effect of the AB potential on the stability of the system qualitatively. From Eq. ([5]) we see 
that the effect of the AB potential appears only in increasing the angular momentum I. Classically it increases the 
centrifugal force on the particle. Thus it tends to counteract the tendency of the particle being pulled towards the 
negative continuum. Hence one expects that the AB potential ([3]) will stabilize the system against pair production. 

Quantitatively, to determine the critical charge a cr {B) in the presence of the AB potential, it is sufficient to consider 
the range of electron energies near the negative-energy continuum limit —to. Introducing functions F(r) = r/(r) and 
G(r) — rg(r), and eliminating G(r) from ([5]), we obtain the equation for the function F(r) with E rj —to in the form 

d 2 F ( 2 2 2Ea a 2 - (I + e B)(l + e B + 1) 

, + IE — to H 1 = 

dr z \ r r A 

The solution for G(r) near E = —to can be found from the equation 

n , . l + eoB + 1 rdF 

G(r) = F(r) — s — — . 

a a dr 

The solution of Eq. (f27|) which tends to zero as r — ► oo can be expressed through the Whittaker function of the 
form 



F = 0. 



(27) 



F(r) ~ W U g{2\r), (29) 

where 

Ea 



/?=— , = v/ a 2 _ (/ + eo B + 1/2) 2 (30) 



near E = —m, or through the MacDonald function of imaginary order 



F(r) - v^WV8mar) (31) 

for E = —to. It follows from Eqs. (|29[) or (|3"Tj) that the bound electron state with E w —to is localized in space. 
Such behavior of the electron state can be easily explained if we treat Eq. (|27f as the Schrodinger equation with the 
effective energy e = (E 2 — m 2 )/2m and the effective potential 

2Ea a 2 - {I + e B)(l + e B + s) 
u eff(r) = r - 2 • (32) 

In the case E w —to for the ground electron state I = the potential U e ft(r) has a form of a wide barrier. One notes 
that at large distances from the Coulomb center for an electron with energy E w — m the effective potential is not 
attracting but repulsive. In the presence of the magnetic flux the height and width of the effective potential barrier 
increase for s = 1 and e$B > 1. Therefore, the probability of pair production by the Coulomb field decreases in the 
presence of the magnetic flux. 

For a weak magnetic flux the form of vector Coulomb potential for r < R is not essential to the principal result. 
The calculation is most easily performed if we consider the simplest model with only the 3- vector Coulomb potential 

e it 
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in the range r < R. Then the radial solution F(r) that is finite at r — in the range r < R is expressed via the Bessel 
function of integer order \l\ as 



F(r)~rJ|,|(cr), (34) 



where 



Applying the continuity relations 



£+-^) 2 -™ 2 . (35) 



G(r)\ (G{r) 



(36) 



and taking into account the fact that the parameter R must be small compared with 1/m and also that E w — m, we 
obtain the transcendental equations (for I = 0) that implicitly determine the energy spectrum of an electron near the 
negative-energy continuum boundary -mas 

cB M = 1 - L^4±) , (37) 

and the critical charge for the ground state 

2MB) «M^ = 1 _(,SM) . (38) 

Ma cr {B)) V K iv [z) J z=y / Smaar(B)R 

Here 



/9 = "X' ^V^M+W, I / = 2 A / a 2 r (B)_( eo B + 1/2)2, (39) 

and the prime denotes differentiation with respect to the argument x of the Whittaker function in (|37[) . or the argument 
z of the MacDonald function in (|38|) . 

Equations (|3"Tj) and (|38[) can only be solved numerically. These equations can be simplified somewhat for Rm <C 1 . 
For instance, to simplify Eq. (|38p . we represent the MacDonald function in the form 

^ = 2^)[ J --^' (40) 



where 



EI /z\c+ 2<; 
fciroT+fc + i) \2) ■ ( ] 



, =o fc!r(/i + fc + 1) 



Keeping only the lowest order terms in the expansion of the MacDonald function at small values of the argument, we 
obtain 

/ \ 1/2 

K iv (z) ~ - , ^ . sin (i/ MM/2) +ar g r(l-ii/)), (42) 
\^smn(i/7r) y 

where T(z) is the Euler gamma function. With Eq. (|4"2")> . we finally obtain an approximate form of Eq. ([55)) at small 
Rm: 

2a crb Jl ) acrb \ = l-i/cot[i/ln(|«|/2) + argr(l-ii/)]. (43) 

Jo{ a crb) 

Here a cr b = a cr (B). 

Numerical solution of Eq. (|43|) gives a cr b ~ 0.79 for B — 0, a cr i, s» 0.84 for eo-B = 0.1, a cr b s* 0.91 for eo-B = 0.2 
at i?m = 0.02, and a crb « 0.70 for 5 = 0, a cr6 « 0.77 for e B = 0.1, a crb w 0.84 for e S = 0.2 at i?m = 0.006. 
One sees that the critical charge increases with the increase of B and decreases with the decrease of Rm. Therefore, 
the instability of the vacuum of quantum electrodynamics in a strong vector Coulomb field in 2+1 dimensions in the 
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presence of the magnetic flux tube of a very small radius must occur at larger critical charge, as compared to that in 
the absence of the magnetic flux. Thus, the magnetic flux stabilizes the vacuum of quantum electrodynamics in the 
vector Coulomb field. 

Similarly, Eq. (|57|) can be simplified. With the condition Rm « 1 we use the representation of the Whittakcr 
function for small values of the argument in the form 

Wp,teW - r(l/2-/3 + *0) + T(l/2-f3-i9) ' [ ' 

= \ T ( 2l9 )\ 2:z; 1/2 cosf$(Vn (45) 

\T(l/2-0 + iff)\ cos m^> I4 j 

where 

=-01nx + argr(2i0)-argr(l/2-^ + i6>). (46) 

Taking the derivative with respect to x in Eq. ([45)1 and substituting the resulting expression and function ()45[) in Eq. 
([37[1 . we finally obtain the simplified transcendental equation 

-6ki(2\R) + avgT(2i9) - ar g r(l/2 - (3 + iff) = tan" 1 Y + im r , (47) 

where 

y = *-i(!-cR#^y (48) 

V2 UcR)) V ' 

The equation for the energy spectrum for any integer I can be derived by the analogous way. 

Equation ([47|) with a < a cr b implicitly determines the eigen spectrum of bound electron states for I = with 
m > E > —to. It can be shown there are real solutions of this equation only for a < a cr bi but with a > a cr b there 
is a formal solution of the form E = E — iw, where Eq = — m — ci(a — a cr b), c\ ~ 1 for the lowest state. With 
a — a cr b <C a cr b, the imaginary part w is exponentially small. Such a solution can be found from formula (j47[) through 
analytic continuation of E as a function of a into the range a > a cr b- However, for a — a cr b -C a cr b, the imaginary 
part can be more readily determined in another way. Indeed, the appearance of the imaginary part means that for 
a > a cr b, the corresponding Dirac equation has only a formal solution with E = Eq — iw for the electron states. 
However, for a > a cr b, the same equation also describes a positron with the energy Eq = to + c\(a — a cr b) since the 
Dirac equation for a positron for a > a cr b can be obtained from the Dirac equation for an electron by replacing E 
with —E and a with —a. Hence, for a > a cr b, the positron states are quasi-stationary. For a — a cr b = Aa <g; a cr b, the 
width of the quasi-stationary level w can be estimated in the semiclassical approximation. To obtain this estimate, 
we must compute the transmission coefficient through the potential barrier in Eq. I|27p. We note that the width 
w is half of the reciprocal of the positron lifetime or twice the probability of pair creation by the Coulomb field. If 
Aa <C a cr b, the positrons created by the field are very slow, and the Coulomb barrier is hardly transparent for them. 
The probability of pair creation is therefore exponentially small, i.e., 



w ~ m exp 



-2na 



Ve 1 



mexp ( -C2\j-^- ) , c 2 ~ 1. (49) 



IV. SUMMARY 



In this paper we solve exactly the Dirac equation in 2+1 dimensions in a combination of Lorentz 3- vector, scalar 
Coulomb and AB potentials. We employ the two-component Dirac equation which contains a new parameter intro- 
duced by Hagen to describe the spin of the spin-1/2 particle. We study the energy spectrum of an electron near the 
negative-energy continuum boundary and the critical charges for some electron states. Fermion pair production from 
a vacuum by a strong Coulomb field in the presence of the magnetic flux tube of zero radius is also considered. It is 
shown that the presence of the AB flux tends to stabilize the system. 

Finally, we note that solutions to the two-component Dirac equation in the AB potential coincide with the solutions 
of the Dirac equation in 2+1 dimensions for a massive neutral fermion with the anomalous magnetic moment in a point 
charge placed at the origin z = (see, for example, Ref.[20]). In the three-dimensional space, such a field corresponds 
to the electric filed of a thin thread that is perpendicular to the plane z = and carries the electric charge with the 
constant linear density. Thus, the solutions to the two-component Dirac equation in the AB potential can be directly 
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applied to the planar scattering of a massive neutral fermion with anomalous magnetic moment interacting in the 
electric field of the thin thread, which was first predicted by Aharonov and Casher in Ref. plj . 
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